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We reanalyze the Fermi spectra of the Geminga and Vela pulsars. We find that 
the spectrum of Geminga above the break is exceptionally well approximated by 
a simple power law without the exponential cut-off, making Geminga's spectrum 
similar to that of Crab. Vela's broadband gamma-ray spectrum is equally well 
hM fit with both the exponential cut-off and the double power law shapes. 

pH In the broadband double power-law fits, for a typical Fermi spectrum of 

Q-i a bright 7-ray pulsar, most of the errors accumulate due to the arbitrary 

O parametrization of the spectral roll-off. In addition, a power law with an ex- 

ponential cut-off gives an acceptable fit for the underlying double power-law 
spectrum for a very broad range of parameters, making such fitting procedures 
insensitive to the underlying Fermi photon spectrum. 
>■ Our results have important implications for the mechanism of pulsar high en- 

ergy emission. A number of observed properties of 7-ray pulsars, i.e., the broken 
00 power law spectra without exponential cut-offs and stretching in case of Crab 

beyond the maximal curvature limit, spectral breaks close to or exceeding the 
(^ maximal breaks due to curvature emission, a Crab patterns of relative intensities 

^ of the leading and trailing pulses repeated in the X-ray and 7-ray regions, all 

IJ point to the inverse Compton origin of the high energy emission from majority 

•^ of pulsars. 

1. VERITAS detection of Crab pulsar: a case for inverse Compton scattering 

origin of 7-ray emission 

The recent launch of the Fermi Gamma-Ray Space Telescope and subsequent detection 



of a large number of pulsars (Abdo 2010d) revolutionized our picture of the non-thermal 
emission from pulsars in the gamma-ray band from 100 MeV up to about 10 GeV. At even 
higher energies, in the very-high energy (VHE) band, the detection of the Crab pulsar at 



25 GeV by the Magic Collaboration (Aliu & MAGIC Collaboration 2008) and recently at 



120 GeV by the VERITAS Collaboration (IVERITAS Collaboration & Aliu|[20ll| in the 



very-high energy (VHE) band allow to stringently constrain the very-high-energy emission 
mechanisms in the case of the Crab pulsar. 



Lyutikov et al. (2011) have argued that in case of Crab the inverse Compton scattering 



is the main emission mechanism of the very high energy emission. This was based on the de- 



tection of Crab pulsar by VERITAS collaboration above 150 GeV (VERITAS Collaboration 



& Aliu 2011) with non-exponential cut-off above the spectral break. The non-exponential 



cut-off was later confirmed by the MAGIC collaboration (Aleksic 2011 which also preferred 
IC model over the curvature emission). 

The curvature emission in pulsars is limited to energies below 
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where Rl is the light cylinder radius, P is pulsar period of rotation, ^ is a dimensionless 
scaling parameter ^ = Rc/Rl, Rc is the radius of curvature of magnetic field lines, B = 
Bns{Rns/R)^, where B^s is the magnetic field on the surface of the neutrons star and Rns 
the stars surface and rj = E/B < 1 is the relative strength of the accelerating electric field , 



(Lyutikov et al. 2011) 



If the 7-ray photons are due to the curvature emission of a radiation reaction-limited 
population of leptons, the spectrum above the break must show an exponential cut-off. The 
detection of the Crab pulsar by VERITAS collaboration ( VERITAS Collaboration fc Aliu] 
2011) clearly demonstrated the non-exponential cut-off above the spectral break, see Fig. [Ij 



Lyutikov et al. (2011) have argued that this is inconsistent with the curvature emission. 



2. Spectral breaks in 7-ray pulsars 

The curvature and the IC model of the high energy emission offer different interpreta- 
tions of the spectral breaks. In case of the curvature emission, the spectral break corresponds 
to the maximal energy of the radiation reaction-limited acceleration. Above the break the 
spectrum is expected to have an exponential cut-off. In the IC model the spectral break 
corresponds to the break in particle distribution function, which is likely to be of the power 
law type both below and above the break. Thus, observations of the non-exponential cut-off 
above the break favor the IC model. 

Crab pulsar is a bright emitter in all spectral bands. Thus, there is a lot of low frequency 
target photons available for Compton scattering. What about other pulsars? In the present 
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Fig. 1. — High energy spectrum of Crab demonstrating the non-exponential spectral break 



inconsistent with curvature emission (figure from VERITAS Collaboration & Aliu 2011) 



paper we argue that that there are evidence for a universal dominance of the inverse Compton 
scattering mechanism in 7-ray pulsars. 



First, Lyutikov et al. (2011) compared the observed spectral breaks of Fermi pulsars 



from the first Fermi catalogue (Abdo 2010d) with the predicted breaks due to curvature 



emission, Eq. ([T| by calculating the ratio of the observed Ehr and predicted spectral break 
ebr, see Fig. [2j A significant number of pulsars the ratio is close to one and for one pulsar, 
PSR J1836 + 5925, the ratio is even larger than one. In order to explain the spectral break 
for these pulsars as a result of curvature radiation, an accelerating electric fields should be 
close to or even larger than the magnetic fields . 

What is more, the example of Crab demonstrates that the spectral break may not be 
related to the maximal curvature photons. The presence of non-exponential breaks in pulsars 
others than Crab would be a clear indication that the break is not related to maximal energy 
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Fig. 2. — Ratio of the observed break energies Eh,, for 46 pulsars to the niaxinium predicted 
for curvature radiation ebr, which is given by Eq. ([I| with r] = ^ = 1 (Lyutikov et ah 2011). 



of curvature emission by radiation-hmited acceleration of leptons. Realistically, due to low 
photon counts, only the brightest 7-ray pulsars - Crab, Vela and Geminga - allow precise 
enough measurements of fluxes to distinguish between different spectral shapes. In addition, 
typically, only phased average properties have enough photon counts to distinguish between 
the models. 



3. Spectra of Geminga and Vela pulsars 
3.1. Geminga: a non-exponential cut-ofF 

Geminga is one of the brightest 7-ray pulsars. Its spectrum is, conventionally, fit with 



a power-law plus exponential cut-off (Abdo 2010a). We have performed an independent fit 



to Geminga spectrum, see Table M and Figs. ^W5] For the fits, the processed data (energy 



flux) provided by the Fermi collaboration were used, Table 3.1 



First, we performed ~x^ fits of the whole spectrum using a particular type of double 
power-law prescription (with 4 independent parameters, row (a) in Table |2| the power-law 
plus exponential cut-off (with 3 independent parameters, row (b) in Table [2]) and a softened 
exponential cut-off (with 4 independent parameters, row (c) in Table [2]), see Table [2] and Fig. 
3J Weighted and unweighted reduced y^ fits were not statistically distinguishable. (Overall 
normalization is the additional fit parameter to the ones listed in Table pi) 
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Fig. 3. — Fits to Geminga spectrum by models described in Table [2j Solid line is double 
power-law model (a) from Table [2] (reduced y^ = 1.26), dot dashed is an exponential cut-off 
(reduced x^ = 1-13) and dashed is a softened exponential cut-off (reduced x^ = 0.83). Since 
in all cases x^ is of the order of unity, all three models are statistically acceptable. (See also 
Fig. lilfor the error analysis of the double power-law fit.) 



Inspection of Table M and Figs. K3M tells, first, that since the reduced x is close to 
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Table 1: Phase- averaged data for Geminga pulsar. 
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Table 2: Fits to the broad band spectrum of Geminga. The values for fit (c) can be compared 



with the best fit done by Fermi team, /3 = 1.12, etr = 1-58, b = 0.81, (Abdo 2010a). 



unity all three models are statistically acceptable. (In Appendix p^ we also demonstrate that 
for a typical Fermi spectrum of a bright 7-ray pulsar a fit with an exponential cut-off can 
actually give a satisfactory approximation to the underlying power-law spectrum for a very 
wide regime of fit parameters, e.g., over two orders of magnitude in break energy.) 

Second, the fit that looks the best down to low fiux levels, fit (a), actually has largest 
reduced x^ among the models. The errors for the this fit are accumulated at intermediate 
energies due to an arbitrarily chosen parametrization of the spectral roll-off. Since we are 
mostly interested in the high energy scaling of the spectrum, we can take only the highest 
energy points. Having chosen the last 6 points we derive the fit given in Table |3j see also 
Fig. i 
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reduced x^ 
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Table 3: Properties of the fit of the high energy tail of Geminga. 



Thus, the high energy part of the Geminga spectrum is exceptionally well fit with a 
power-law function, giving reduced x^ = O-l- This exceptional small value of x^ is a bit 
surprising: despite the lower overall count the highest energy data points have the smallest 
error bars, Table [3TT| due to the very low background counts above ~ 2 GeVs. A small value 



of x^ for the high energy fit may be due to the fact that the fit has only four degrees of 
freedom - in this case the random chance for all the points being very close to the model is 
not small. This may also be an indication that the number of energy channels was chosen 
not in an optimal way, e.g., too broad energy bins. (If only 5 highest energy points are used, 
the fit values are a = 2.8, x^ = 0.03 for 3 degrees of freedom.) 
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Fig. 4. — Errors for the broadband fit to Geminga data using double power law function. 
The lower panel shows the residuals (on a hnear scale) for the double power-law fit, while the 
upper panel shows relative errors that are used to calculate x^- The errors are not random 
indicating the inadequacy of the arbitrarily chosen spectral roll-off. Also, most of the x^ 
is accumulated near the break energy due to the arbitrary parametrization of the spectral 
roll-off. Note that the highest energy data points actually have the smallest error bars. 

3.1.1. Geminga to be seen at very high gamma-ray energies? 
The spectral index above the break dN/dE oc -E~^ with p fa 5 for Geminga is steeper 



than that of the Crab pulsar [p = 3.8 above the break McCann & VERITAS Collaboration 



2011 ). Geminga is nearly two and a half times brighter than Crab at the peak of its spectral 
distribution of ~ IGeV, but its has lower break energy and most importantly the steeper 
spectrum. The expected flux at ~ 120 GeV, the low energy threshold for the VERITAS 
Cherenkov telescope, is then few ~ 10^^ MeV cm~^, s~^. This is approximately ten times 
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5. — Fits to the high energy tail of the Geininga spectrum: power law (solid line, 
0.1) and exponential cut-off (dashed line, x^ = 2). 



lower than the flux from the Crab pulsar. Geminga, still, has a chance to be detected, since 
unlike Crab it has no strong background from the PWN. (It is not clear whether the excess 



TeV flux from the general direction of Geminga (Abdo et al. 2008) is related to the pulsar 



see, though Salvati & Sacco 2008)). 



3.1.2. Geminga: Phase-resolved fit 

Phase variation in the cut-off energy can affect the overall spectral fit, but not at the 
highest energies. For example, Abdo (2010a) cite the cut-off energy variations between 



approximately 1 and 3 GeV, while our high energy fit starts at 8.5 GeV, nearly three times 
higher, see Table |3] and Fig. [5j In addition, due to a low photon count in the phase- resolved 
SEDs, they naturally miss the highest energy, lowest counts tails. 

To test the spectral shape in phased-resolved spectra we have performed spectral fitting 
at one particular phase, 0.614 < (j) < 0.623, corresponding, approximately to the second 
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peak Abdo (2010a). Since phase resolved data were not available, we assumed the error on 
the flux to be equal to 2 x 10~^^ergcm~^ s~^, see Table E^ and Fig. m 
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Table 4: Parameters for the phase-resolved fit of Geminga. 

Thus, both the double power-law and the exponential cut-off provide a statistically 
viable description of data for a particular phase-resolved spectrum of Geminga. 



3.2. Vela: consistent with the non-exponential tail 

Vela, the brightest 7-ray pulsar, offers tantalizing, but non-conclusive evidence of the 
non-exponential break, see Fig. [7| We performed an unweighted x^ fit with phase-averaged 
data provided by the Fermi collaboration. Table |3.2[ We used various double power-law 



prescriptions (with 4 independent parameters), the power-law plus exponential cut-off (with 
3 independent parameters) and a softened exponential cut-off (with 4 independent param- 
eters), see Table |6} Spectrum of Vela is equally well fit with a softened exponential cut-off 
and a double power-law spectrum. 

In case of Vela, a particular choice of the parametrization of the double power-law 
spectrum illustrates that an arbitrary choice of the roll-off between the two spectral power- 
law components has an important effect on the fit: while the more conventional double 
power-law parametrization in line (a). Table [61 gives unacceptable fit, double power-law 
spectrum parametrization in line (b) in Table [6] gives an acceptable fit. 



4. Optical- V-ray-7-ray correlation 

Within the framework of the SSC model the power emitted by IC is related to the power 
of the seed photons. Photons of different energies that are emitted by the same particles 
should in principle produce similar pulse profiles. In our model one expects, therefore, that 
the pulse profiles in X-ray and in gamma-rays are similar because the secondary plasma 
emits synchrotron radiation in X-rays and IC scatters UV photons into the VHE band. And 
indeed, the ratio of the amphtudes of the two pulses in the pulse profile of the Crab pulsar 
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Table 5: Phase-averaged data for Vela pulsar. 
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Fig. 6. — Phase-resolved fit to Geminga spectrum using double power law (solid line, reduced 
X^ = 0.46) and exponential cut-off (dot dashed line, reduced x^ = 0.23). The fit corresponds 
to the phase 0.614 < < 0.623 of lAbdol (]2010al). 
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Table 6: Parameters of the spectral fit for Vela pulsar 



changes consistently in the X-rays / soft gamma-ray band and in the high energy gamma- 
ray band, see Fig. |8} In X-rays the main pulse dominates over the inter pulse. The ratio 
changes towards higher energies and reverses in the soft gamma-ray band at about 1 MeV. 
Mirroring the evolution at lower energy, the main pulse dominates at 100 MeV, while at 
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Fig. 7. — Fits to Vela spectrum by models described in Table [6j Solid line: model (a), 
dotted line: model (b), dashed line: model (c), dot dashed line: model (d) from Table [6] 
Models (a) and (b) are double power-laws, (c) is an exponential cut-off and (d) is a softened 
exponential cut-off. Models (a) and (c), as well as models (b) and (d) have approximately 
the same reduced y^ . 



120 GeV the interpulse clearly again dominates over the main pulse (Abdo||2010d VERITAS 



Collaboration & Aliu 2011). 



Vela also offers tantalizing evidence for the correlation between the optical/X-ray and 
7-ray bands, Fig. |9| all three 7-ray peaks (PI, P2 and P3 in the top panel) have lower 
frequency analogues: PI in the X-rays, P2 in optical and P3 in both bands. This is expected 
in the IC model, though relative intensities and spectra of the components depend on the 
details of the photon and particle distribution. 

Crab pulsar is somewhat exceptional among the brightest X-ray-7-ray pulsars since it 
has no indication of the thermal component. The X-ray emission of Vela and Geminga 



are strongly dominated by the thermal components coming form the surface (Halpern & 



Ruderman 1993 Pavlov et al. 2001 Manzali et al. 2007). In Vela, the non-thermal component 
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Fig. 8. — Evolution of the Crab profile with energy ( Abdo|[2bl0c ). Note that that the lower- 
energy evolution of the increasing interpulse to main pulse ratio is mirrored in the 7-rays. 



Such behavior is expected in synchrotron-self-Compton model (Lyutikov et al. 2011). 



in X-rays is barely detected at some orbital phases, roughly aligned with the 7-ray peaks, 



especially with P2, see Fig. 9 of Manzali et al. (2007). The relative weakness of the non 



thermal X-ray components does not allow a measurement of the evolution of the P1/P2 
ratio with energy, like in Crab. Overall, both Vela and Geminga spectra show a complicated 
phase- and energy-dependent mix of thermal and non-thermal components {e.g., Kargaltsev 



et al. 2005 Romani et al. 2005). In both cases the optical emission, the main target for IC 



scattering, is highly non-thermal. 
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Fig. 9. — Vela profile in broad energy bands (Abdo 2010b). Chandra profile in the soft 



X-rays is dominated by thermal emission and is not of interest here. Optical and RXTE 
signals are mostly non-thermal. Note that all three 7-ray peaks (PI, P2 and P3 in the top 
panel) have lower frequency non-thermal analogues: PI in the X-rays, P3 in optical and P2 
in both bands. This is expected in the IC model, though the spectra and relative intensities 
depend on the details of the photon and particle spectra. 

5. Discussion 



In this paper we demonstrate that, first, broadband pulsar high energy spectra are gen- 
erally consistent with the double power law shapes. This results is based on three brightest 
7-ray pulsars. In Geminga the high energy part of the SED is exceptionally well fit with a 
power law, while the spectrum of Vela is inconclusive. Importantly, in Geminga most errors 
in the broadband double power law fit are accumulated due to an arbitrary parametrization 
of the spectral roll-off between the two asymptotic power laws. Thus, the standard x^ fit 
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underestimates the goodness of fit in this particular case. In addition, the low background 
at high energies, above 2 GeV, make Fermi data especially sensitive (smallest error bars) to 
the shape of the high energy spectrum. 

If pulsar spectra generally show non-exponential cut-offs, a fact well established in the 



Crab pulsar by VERITAS (VERITAS Collaboration & Aliu 2011) and in Geminga in the 



present paper, this has important implications for models of pulsar 7-ray emission. This 
implies the importance of the inverse Compton (IC) scattering for the production of 7-ray 



photons (Lyutikov et al. 2011) and would signify a paradigm shift in the study of pulsar 



high-energy emission. 



Lyutikov et al. (2011) argued that in Crab and in a number of other pulsars the inverse 



Compton scattering is the main emission mechanism of the very high energy emission. The 
authors tried to give the most general arguments in favor of IC, mostly independent of 



the numerous possible particular details of a full radiative model. Overall, Lyutikov et al. 



(2011) adopted the Cheng et al. (1986) paradigm of pulsar high energy emission, but with 



inverse Compton scattering playing a more important role that previously assumed. Recall, 



that Cheng et al. (1986) argued that IC scattering is not dominant in Crab, but may be 



important for Vela. VERITAS and MAGIC results on Crab indicate that IC is the dominant 
mechanism in Crab. 



In the model of Lyutikov et al. (2011) the inverse Compton scattering occurring in the 



Klein-Nishina regime by the secondary particles results in a picture that, overall, is consistent 
with observations without any fine-tuning. The key features of our model are (i) A population 
of primaries that is accelerated in a modest electric field, which is a fraction 77 of the magnetic 
field strength near the light cylinder with a typical value of rj is 10~^. The suppression of the 
scattering cross-section in the Klein-Nishina regime (and the corresponding lower radiation 
loss rate of electrons) allows primary leptons to be accelerated to very high energies with 
hard spectra, (ii) The gain in energy of the primaries in the electric field is balanced by 
similar curvature radiation and IC losses (radiation reaction limited); (iii) The secondary 
plasma is less energetic, but more dense and has approximately the same energy content as 
the primary beam. The secondary plasma is responsible for the soft UV-X-ray emission via 
synchrotron/cyclotron emission and the high energy 7-ray emission that extends to hundreds 
of GeV via the inverse Compton process. The IC emission from the primary beam extends 
well into the TeV regime but will be difficult to detected due to the low predicted fluxes. 

The spectra of all three brightest 7-ray pulsars are consistent with the broken power- 
law distribution (and, thus, with the inverse Compton origin of the 7-ray emission). What 
is more. Crab and Geminga are inconsistent with exponential cut-off (and, thus, with the 
curvature origin of the 7-ray emission). Thus, inverse Compton scattering may be the 
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dominant radiation mechanism in all pulsars, though a curvature contribution in some pulsars 
cannot be excluded. 

In the analysis we use the processed data points provided by the Fermi collaboration as 
well as data from the published papers. These data points are actually somewhat dependent 
on the model used to derive them, since the model of the spectrum goes into untangling 
the data from the instrument response. This introduces some uncertainty, but since we are 
looking at the brightest pulsars we expect that a particular model used to estimate the flux 
in individual spectral bins is not going to make too much of a difference. We would like to 
encourage the community to test our results using the raw count rates data. 

It is not clear what parameters determine the dominance of the IC scattering among 



Vela, Crab and Geminga. Calculations of Cheng et al. (1986) suggested that IC scattering 



may be important for Vela, but our results suggest that Vela is the weakest case for the IC 
among the three. The case for dominance of IC scattering in Crab was, partially, expected 
due to high soft photon luminosity (e.g., Hirotani]|2007). What is more, the Crab emission is 



completely non-thermal and originates in the same region of the magnetosphere as the 7-ray 
emission, while in Vela and Geminga a large contribution to X-ray emission comes from the 
thermal emission from the surface. 

In addition, the strong energy dependence of IC scattering in the KN regime favors the 
lower energy optical-UV photons. In this view, the example of non-exponential cut-off in 
Geminga, and the implied importance of the IC scattering, are quite surprising given the low 
observed optical-X-ray luminosity. This poses challenges for radiative modeling. It is likely 
to be related to the fact that inverse Compton scattering is highly dependent on the details 
of both the soft photon distribution and the particle's distribution {e.g., their anisotropy). 

Finally, let us recapitulate the main arguments in favor of inverse Compton scattering 
being the dominant pulsar radiation mechanism at GeV energies: 

• The spectrum of Crab pulsar extends beyond the upper limit for curvature emission 

• The beak energy in many pulsars approaches and even exceeds the upper limit for 
curvature emission 

• In the three brightest pulsars. Crab, Geminga and possibly Vela the spectrum above the 
break is not exponentially suppressed, implying that the break is not due to curvature 
remission of leptons in the radiation reaction-limited regime 

• In Crab, the energy dependence of the relative intensity of emission peaks at 7-ray en- 
ergies mirrors the lower energy energy dependence, as expected in the Inverse Compton 
model. 



I would like to thank Oleg Kargaltsev, Matthew Lister, Rafael Lang, George Pavlov, 
Scott Ransom and Mallory Roberts for help with the data analysis, Roger Blandford, John 
Finley, Andrew McCann, Roger Romani, Nepomuk Otte, Mark Thieling and David Thomp- 
son for comments on the manuscript and Fermi Large Area Telescope Collaboration for 
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A. On the exponential fit to Fermi-like data with power-law distribution 



Often, Fermi pulsar data are fit with a softened exponential function, oc e"e' 



-(f/^tr) 



which is supposed to imitate variations in the exponential cut-off energy. In this appendix 
we address a question: if the data come from double power-law distribution, how well will 
they be represented by the softened exponential fit? As an example, we fit a function 
/ = 2/(x + 1/x) sampled at equally spaced logarithmic intervals from 0.1 to 40 with a 
typical error of 0.01 and random scatter 0.01. This choice of the parameters simulates a 
typical bright pulsar spectrum sampled in unites of GeV and maximum flux of the order of 
~ 100 of the statistical error. 

We find that the fit function provides an acceptable fit (reduced x^ ^ 1) ^oi an excep- 
tionally wide range of the break energies 0.01 < ebr < 2, see Fig. [lOJ For less bright pulsars, 
where due to low photon statistics the spectrum can be measured up to a lower maximal 
energy, the reduced x^ remains smaller than unity even for larger values of the fitted ef,r ■ 

Thus, if a typical spectrum of a bright Fermi pulsar is a broken power-law, the spectral 
fit with a softened exponential function, oc e^e"^'^/'^'"'-' , gives acceptable fits for a very wide 
range of parameters and, thus, does not provide a sufficiently sensitive probe of the underlying 
spectrum. 



This preprint was prepared with the AAS I^Tj^X macros v5.2. 
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Fig. 10. — Numerical fit to the function / = 2/[x + 1/a;) using a softened exponential 
function, oc e"e~^'^''^*'^'' . Two curves correspond to extreme values of the parameters {e^r = 
0.01, a = 1.87, b = 0.35 and thr = 2, a = 1.68, b = 0.98), but both give acceptable fits 
(reduced x^ = 0-33 and 1.01 correspondingly.) 



